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. Abstract. The correlation often observed in blazars between optical-to-radio outbursts and 

pl"^ ' gamma-ray flares suggests that the high-energy emission region shall be co-spatial with the 

• , radio knots, several parsecs away from the central engine. This would prevent the important 

i-^ ' contribution at high-energies from the Compton scattering of seed photons from the accretion 

disk and the broad-line region that is generally used to model the spectral energy distribution of 
^ , low-frequency peaking blazars. While a pure synchrotron self-Compton model has so far failed 

i-^ • to explain the observed gamma-ray emission of a fiat spectrum radio quasar like 3C 279, the 

I inclusion of the effect of multiple inverse-Compton scattering might solve the apparent paradox. 

■ Here, we present for the first time a physical, self-consistent SSC modeling of a series of shock- 

waves in the jet of 3C 279. We show that the analytic description of the high-energy emission 
from multiple inverse-Compton scatterings in the Klein-Nishina limit can fairly well account for 
^ ■ the observed gamma-ray spectrum of 3C 279 in flaring states. 
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Q^ ■ 1. Introduction 

i The radio-to-infrared emission of blazars is well understood as being synchrotron emission 

associated to bright structures observed to propagate — often with apparent superluminal 
motion — in a relativistic jet at parsec-scale distances [H- The preferred model for the emission 
of these knots are shock waves resulting from disturbances at the base of the jet that move 
^ ■ outward, become supersonic, compress the gas locally, and accelerate particles at the shock 

front [2]. In order to test the validity of this shock-in-jet model, a tool was developed to fit the 
multi-frequency light-curves of the bright quasar 3C 273 with a series of model outbursts 
The same approach was then also used to describe the flaring behavior of other sources: the 
blazar 3C 279 [5], and the micro-quasars GRS 1915-M05 [g] and Cyg X-3 0; [l]. While these 
studies still relied on model parameters that are observables, we now developed a fully physical 
parametrization of the synchrotron emission of shock waves [oj . 

The correlation often observed in blazars between optical-to-radio outbursts and gamma-ray 
flares suggests that the high-energy emission is closely linked to the synchrotron component 
(e.g. [lC|). If the emission arises parsecs away from the central engine, the seed photons for 
inverse-Compton scattering cannot be of external origin. However, for a low-energy peaking 
blazar like 3C 279, it has been shown that such a synchrotron self-Compton (SSC) scenario is 
not able to account for the observed gamma-ray emission [TH : 12]. While this is true for single 



scatterings, the effect of multiple-inverse Compton (MIC) orders has been reconsidered recently 



to solve this apparent paradox |l3l | . We present here a first attempt to apply an analytical model 
for multiple synchrotron self-Compton (MSSC) emission of shock waves to a broad dataset of 
3C 279. 



2. Theory 

This section aims at linking the theoretical model for MIC scattering 13j; [ij] with the parameters 
and the formalism used for the synchrotron emission of shock waves in a relativistic jet [9|. As 
a starting point we link the optical depth tq of ^] with the normalization K of the electron 
energy distribution, assumed to be of the form Uei'j) = K 7"^ for 7min < 7 < 7maxi where ng is the 
electron number density and 7 is the Lorentz factor of the relativistic electrons. For p > 2 and 
7max ^ Trnin, integrating ne(7) leads to K = (p — l)7^7n "-e such that Eq. (7) of becomes: 

where ctt is the Thomson cross-section and xr = JrR is the thickness of the emitting slab 
containing the shocked plasma behind the shock front, which is assumed to be a constant 
fraction fji of the jet radius R 

For a given electron energy distribution ne(7), the MIC spectrum is determined by two main 
parameters a and 7cooi- The spectral index a [Fj^ oc connects the u maxima of the 

synchrotron and the inverse-Compton components. The Lorentz factor 7cool corresponds to the 
energy of electrons that have just enough time to cool down while crossing the emission region, 
i.e. x(7cooi) = /5reicicooi(7cooi) = xr, where /3reic is the average speed of electrons along the jet 
relative to the shock front. In the MIC scenario, there are three distinct cases 



(2) 
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In cases A and B, a = (p — l)/2 coincides with tq = 1 and we note the corresponding value of 
Tcool as 7cooi 0- can derive a simple equation linking 7cooi to a as: 
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where y = mc^ / (hus) and ub = eB /{2iTmc) is the cyclotron frequency for an electron mass, m, 
and charge, e, in a magnetic field, B. Eq. ([5|) is valid in case A and also in case C, if in the 
latter the label "0" denotes a = {p — l)/2 rather than tq = 1. For case B, there is an additional 
term (Tmin/Tabs)^""*^"^" on the left-hand side (LHS). This is due to the fact that in this case 
the Compton-scattered flux is dominated by the lowest energy synchrotron photons, i.e. those 
radiating at the synchrotron self-absorption frequency fabs = '^BTabs- 

In Eq. ([5]) we implicitly assumed that 7abs < Tcool- When 7abs > Tcool, the corresponding 
equation is: 
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where Tabs.o is the value of 7abs when tq = 1 and the equation reduces to Eq. ([5]) when 
Tabs ~^ Tcool- For case B, a modification is again needed to the LHS of Eq. ([6]): the ratio 



7cooi/7abs has to be replaced by 7min/7abs- From Eqs. (13) and (18) of [l3(], one notes that the 
right-hand side (RHS) of both Eqs. ([5]) and ([6|) is given by: 

TjuQ 7m7n 2/3rel mc^ (i^eX ^ h 6 ^j-el , . 

y To Ub aTyxjiB^ mca-rXRH 

where Ue = K mc^ j'^^^ / {p — 2) and Ub = B^/{8it) are the energy densities of the electrons 
and the magnetic field, respectively, and where we used Eq. ([1]) for tq. This equation is to be 
calculated for a given jet model. We now have to obtain equations for a that are independent 
of 7cooi in each case A, B and C, and both for 1) 7abs < 7cooi with Eq. ([5]) and for 2) 7abs > 7cooi 
with Eq. ([6|). We get the following equations for a: 

Al) [21ny]a2 + [lnTo + 21n(RHS)-(p-l)lny]a-(p-l)ln(RHS) = (8) 

A2) [2 In (y/7abs)] a^ + [{p-l) In (^lUy) + 2 In (RHS)] a - (p- 1) In (RHS 7LV^o) = (9) 
B) Q = (p-l)/2-lnro/(21n7^in) (10) 

CI) a = ln(7^,V(RHSro))/ln(y7^in) (11) 

C2) [In (2//7abs)] a2 + [{p-1) In (7min7abs) + hi (RHS)] a - (p-1) In (7^n/^o) = (12) 

Pratically, as 7abs is not know a priori, we first derive a for the cases Al, B and CI, which are 
independent of 7abs- We then calculate 7cooi for each of the three cases Al, Bl and CI using 
Eq. [5] and then 7abs given by: 

»+4+a e'^g^jp) xrK 
^^^^ =16^^^^-" ^^^^ 



where we used Eq. (7) of [^] and Eq. (20) of [13i] and where ^^(p) is defined in Eq. (4) of [9[. In 
the special case Bl, we actually first combine Eq. ([5]) and Eq. (fT3]l to get an expression of 7cooi 
independent of 7abs- We can now discriminate between the various cases. If 7abs > 7cooi occurs 
we calculate a for the relevant case (A2, B2 or C2) and then update the values of 7cooi with 
Eq. dS]) and of 7abs with Eq. (fTSl) . 

We note that the three different cases A, B and C concern the inverse-Compton spectrum 
and are distinct from the three-stage evolution of the synchrotron self-absorption turnover in 
the shock model of Marscher &: Gear [2| • The latter is governed by the width x of the emission 
region that can be smaller than xr if synchrotron or inverse-Compton cooling timescales, tcooh 
limit the effective emission region to Xabs = /3reicicooi(7abs) < xr. In the MSSC scenario, 
a^abs = 2;_R,(7cooi/7abs)" and this inverse-Compton cooling dominates over first-order Compton 
or synchrotron coohng [11]. When a > 1, Xabs = a;i?(7nooi/7abs), where 7cooi and 7abs are now 
controlled hy Ur alone (i.e. with Us = in Eq. (8) of [9|). This synchrotron-cooling stage only 
sets in if 7cooi < 7abs is still satisfied when a = 1. Otherwise, we have a direct transition from 
the inverse-Compton stage to the adiabatic expansion cooling stage, where radiative cooling is 
negligible (i.e. Xabs > xr) and thus the emission region is limited by xr. 



3. Method 

In order to test the MSSC scenario for a shock wave propagating in a relativistic jet flow we 
had to fully describe the expected evolution with time of both the synchrotron and the inverse- 
Compton spectra. For this we basically took the same assumptions as described in Sect. 3 of [§] 
to define the parameters relevant for the shock-in-jet modeling. We ended-up with an evolution 
of several key parameters with distance, X, of the shock-wave from the apex of the jet, which is 
itself related to observed time. Those parameters were the normalization, K, and the two limits, 
7min and 7max5 of the electron energy distribution, as well as the strength of the magnetic field. 



B, and the maximal thickness of the emission region, rc/j, which is directly proportional to the 
jet radius R (cf. Sect. [2]). We then calculated the evolution of tq with Eq. ([1]) and derived also 
that of Q, 7cool and 7abs) which are all smoothly changing with time across the transitions from 
one case to the other. 

We then calculated the evolution of the self-absorption flux density Fabs = Fu^^'^^i^ahs) with 
Eq. (6) of 0], which normalizes both the synchrotron and inverse-Compton spectra. The exact 
shape of these components is controlled by p and a for what concerns the spectral slopes and 
by 7mirn 7cooi> 7abs for the position of the various breaks. Depending on the relative order of 
these characteristic Lorentz factors, the synchrotron spectrum can take five different shapes (l5l |. 
whereas the inverse-Compton spectrum has basically three different shapes for cases A, B and 
C [13I ]. The analytical formulation used for the breaks has been simplified as described in 0]. 
A high-energy cut-off is applied to both spectra at the frequency corresponding to 7max- The 
inverse-Compton spectrum is also limited at low-frequency by a break towards a spectral index 
(Fi, oc i''^) of s = {p— l)/2 below the rest-frame frequency i^b^cooI ^^'^ another one to an index 
of s = — 1 below i^B^min' '^■B7abs 7abs > 7min- Finally, we attempted for a smooth transition 
from a multiple-Compton to a first-order Compton spectrum. 
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Figure 1. Fit of ten lightcurves of 3C 279 in 
the optical-to-radio domain with a series of self- 
similar synchrotron model outbursts following 
the spectral evolution shown in Fig. [2j The 
contribution of individual outbursts is shown 
with different colors and vertical lines at the 
top show the times of the available SEDs. 



Figure 2. Spectral evolution of the average 
synchrotron outburst in 3C 279 as derived from 
the fit of 19 lightcurves (vertical lines) as shown 
in Fig. [TJ Synchrotron spectra at different 
times (grey curves) rise and decay with a 
maximum following the red line. The peak of 
individual outbursts is shown by colored dots. 



4. Results 

To test whether the model outlined above for the high-energy emission of blazars is valid for a flat 
spectrum radio quasar (FSRQ) like 3C 279, we took all the spectral energy distributions (SEDs) 
of this object obtained during the operations of the Compton Gamma- Ray Observatory [lit]. 




Figure 3. Spectral decomposition of two SEDs of 3C 279 derived by fitting them self- 
consistently together with ah the optical-to-radio lightcurves (see Fig. [1]). The emission of 
different model outbursts - corresponding to a succession of different shock waves in the jet - 
are shown with the same colors as in Figs. [1] and El The important GeV variability between 
epoch P5b (right panel) and epoch P5a (left panel) in only '^10 days is not fully accounted 
for by the early emission of a rising outburst (green curves), where MIC scattering dominates 
strongly over first-order Compton (dashed line). The points from the ten other SEDs and the 
MAGIC data of 2006 and 2007 [IQ] are shown in grey for comparison. 

We then fitted these eleven SEDs together with 19 lightcurves in the radio-to-infrared domaiifl 
with a set of 24 model outbursts corresponding to a succession of shock waves propagating in the 
jet. Individual outbursts only differ in shock compression factor, rj, and distance from the apex 
of the jet, Xp, needed to build-up the shock [9,]. The obtained fit to the lightcurves is shown 
in Fig. [U while the corresponding synchrotron spectral evolution is shown in Fig. [2j We note 
that the latter has a similar shape to what was obtained previously with only about half of the 
current dataset [![ , but here the rising and slowly decaying peak of the evolution is controlled by 
a relatively high value of 7min and with no synchrotron stage. The overall spectral decomposition 
is illustrated in Fig. [3] with the most extreme high-energy spectrum in the dataset (epoch P5b) 
and the one taken just ~10 days before. In contrary to an early attempt fl^, a pure SSC model 
seems now to be able to account both for the X-ray and high-energy fiaring emission of 3C 279. 
We note, however, that the presented model is not a true one-zone model. Indeed, to have a 
deep enough gap between the synchrotron and the inverse-Compton humps, the latter had to 
be scaled down by a factor 0.1 that corresponds to assuming that the high-energy emission site 
is 10 times smaller than the synchrotron emission region [13]. 

The model we used here assumes a conical jet with a constant Doppler factor, an opening 
radius of 2° pointing 3° away from the line of sight and with a bulk Lorentz factor for the 
shocks of r = 6.8. The main parameters for an average outburst at synchrotron peak fiux are: 
a distance Xp = 16.7 pc from the apex of the jet; a magnetic field of B = 48 mC; an electron 
energy distribution characterized by an index p = 2.10, a normalization Kp = 6460 cm~'^, and 
cuts at 7min = 254 and 7max = 7450. This leads to an excess by a factor of ~ 100 of the 



^ The description of the dataset is postponed to another paper in preparation. 



electron energy density Ue compared to Ub- Although these values are rapidly changing further 
upstream or downstream, they seem reasonable compared to other studies (e.g. [11]). The 
important difference is that our emission site is far more distant from the central engine, where 
the contribution from photons external to the jet becomes negligible. This is at least true for the 
long-term X- and gamma-ray emission. The very early evolution of the outburst in our model 
(see the green spectrum in Fig. [3]) seems to produce flaring emission on time-scales of days at 
GeV energies and up to to the TeV domain as detected by MAGIC on a few occasions [la]. 

5. Conclusion 

We presented a new approach to model the variable high-energy emission of blazars and 
demonstrated its viability by applying it to a very rich dataset of 3C 279 spanning almost 20 years 
and the full SED from the radio to VHE gamma-rays. These preliminary results suggest that a 
shock- in-jet model with a proper treatment of MIC scattering can account for the high-energy 
emission observed in the source at some epochs, while simultaneously reproducing the frequency- 
dependent variability monitored in the radio-to-optical range. Such a pure-SSC scenario for 
FSRQs shall be further tested in a forthcoming paper with the inclusion of the underlying, 
steady jet emission that shall not exceed the observed quiescent state of the source and by 
including addition observational constraints in the X-ray to gamma-ray range, in particular with 
lightcurves of 3C 279 from the satellites Fermi, Swift, and the Rossi X-ray Timing Explorer (e.g. 
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